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ABSTRACT 


Using the quantum Cramer-Rao bound from quantum estimation theory, we derive a funda¬ 
mental quantum limit on the sensitivity of a temperature measurement of a thermal astronomical 
source. This limit is expressed in terms of the source temperature Tg, input spectral bandwidth 
Av, and measurement duration T, subject to a long measurement time assumption TAv » 1. It 
is valid for any measurement procedure that yields an unbiased estimate of the source tempera¬ 
ture. The limit agrees with the sensitivity of direct detection or photon counting, and also with 
that of the ideal radiometer in the regime kTgjhv^f » 1 for which the Rayleigh-Jeans approxi¬ 
mation is valid, where vq is the center frequency at which the radiometer operates. While valid 
across the electromagnetic spectrum, the limit is especially relevant for radio astronomy in this 
regime, since it implies that no ingenious design or technological improvement can beat an ideal 
radiometer for temperature measurement. In this connection, our result refutes the recent claim 


of a r adio astronomy technique with much-improved sensitivity over the radiometer (jLieu et al 


2015D- 


Subject headings: instrumentation: detectors - radiation mechanisms: thermal - methods: ana¬ 
lytical 


1. Introduction 


Astronomical observations are made over frequencies ranging from radio and microwave frequencies 
through to infrared, optical. X-ray, and gamma-ray frequencies. However, measurement techniques and 
instruments vary widely over the electromagnetic spectru m. At radio to millimeter-wave frequencies, the 
chief measurement device in astronomy is the radiometer (Dicke IQdtih . For thermal sources with photon 
flux spectral density no (in photons • s“^ • Hz“^), related to the source temperature Tg, and bandwidth Aiz 
Hz, the sensitivity of an ideal radiometric measurement of no follows the radiometer equation 
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where the left-hand side is the relative sensitivity (Variance / squared mean) of th e radiometer estim ate 

of no and T is the measurement time ( Burke and Graham-Smithiboioll . Recently, Lieu et al. (2015 ) proposed _ 

two -detector setup similar to the intensity interferometer of Hanbury Brown and Twiss ( Hanburv Brown and TwissI 


1957 ) and an estimator of the source temperature that was claimed to achieve the relative sensitivity 
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This equation is essentially the square of eqn. (33) of ( Lieu et al. 2015[) . with the following relabeling for 
consistency with the rest of this paper. We have renamed T of eq. (33) of dlieu et al.lbolill to Tsafyip the 
time needed to obtain one sample in their scheme. Their number of samples N is then N = TjTsamp, where 
here - and throughout this paper - T stands for the total observation time needed to produce all samples. 
We have also replaced their coherence time r with the inverse of the bandwidth Av, which correspondence 
holds to within a constant numerical factor for most thermal spectra of interest. The sensitivity of eq. @ is 
claimed under the condition Tsamp « 1/A^, requiring fast, but not unfeasibly fast, detection electronics. In 
principle however, Tsamp can be made arbitrarily small relative to IjAv « T, so that the second term is the 
limiting one fundamentally. Since hqAvT is the average number of photons incident during the observation 
period, eqn. © shows a Poisson scaling of the relative sensitivity characteristic of coherent-state sources. 
In particular, the sensitivity decreases (i.e., the measurement improves) with increasing source temperature, 
whereas the radiometer equation 0 exhibits no dependence of the sensitivity on uq. In the photon-rich RF 
regime, the scaling of eq. if correct, can offer orders-of-magnitude improvements in sensitivity. 

In (Zmuidzinas 20I5ll. Zmuidzinas has examined in detail the derivation of eo. @, and concluded that 
the two-detector scheme of Lieu et al. ( 20I5ll actually leads to a relative sensitivity 
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(cf. eqns. (136) and (E2I) of (|Zmuidzinasll20I5ri l which is in close agreement with the radiometer equation. 
Moreover, Zmuidzinas has shown that simple photon counting (direct detection) - to w hich the radiome¬ 
ter provides a clos e approximation at radio frequencies for photon-rich thermal sou rces ( Nityanandd llQQ^ : 


ZmuidzinasIl2003a^ - gives an identical sensitivity. In other words, the scheme of ( Lieu et al.ll201^ is no 
bet ter than usual rad iometry and does not provide the purported increased sensitivity. We refer the reader 
to (|Zmuidzinasll2015ll for extensive discussion on the logical fallacy leading to the incorrect result eq. ([2]) for 
the two-detector scheme, and the physical reasons why that scheme cannot work as claimed. 

Our focus in this paper is different. Insofar as the radiation incident on a telescope of any kind, 
irrespective of the wavelength region, is electromagnetic in nature, it is governed by the quantum theory of 
radiation. Any measurement scheme aimed at extracting information about a parameter such as no from 
the incident field is subject to fundamental limitations following from the laws of quantum mechanics. The 
study of the limitations on extracting information fr om quantum systems using quantum measurements 


is the subject matter of quantum es timation theory ( I elstrom 1976 : Holevol 2011 1. which falls under the 


general rubri c of quantum metrology (Giovannetti et al 


information ( Nielsen and Chuana 2000ll . 


201111 with close connections to the field of quantum 


Using the ideas of quantum estimation theory, in particular the quantum Cramer-Rao bound, we show 
that any method of estimating the noise temperature of a thermal source suffers from a minimum relative 
sensitivity equal to the right-hand side of eq. ([3]) in the limit of large observation times T relative to the 
coherence time Tr = IIAu o f the incident radiation. This result provides another refutation, independent 
of the work of ZmuidzinasI (20151. of the claim of eo. ©• More importantly, however, it shows that the 
relative uncertainty of eq. ([3|) cannot be beaten by any conceivable measurement scheme, provided only that 
the scheme yields an unbiased estimate in the sense that the average of many estim ates of the para meter 
converges to the true value - we note that both photon counting and the scheme of (|Lieu et al.ll2015ll lead 
to unbiased estimates and therefore fall under the purview of the bound. 

We briefly me ntion some recent work on quantu m limits for temperature measurement in different 
contexts from ours. iMarzolino and BraunI (|2013l . 12015^ have studied limits on the accuracy of temperature 
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measurement of quantum gases us i ng the quantum Cramer-Rao bound and found that energy measurements 
are optimal. iJarzvna and Zwiera (|2014r) have studied the local estimation accuracy - i.e., the accuracy of 
estimating small deviations of temperature from a preset fiducial value - of pyrometers, i.e., devices that 
measure the total energy output of a thermal source. Under this assumption, they find that the accuracy 
of local estimation of temperature matches that from the quantum Cramer-Rao bound. In this paper, the 
temperature range of the source is not restricted and explicit account is also taken of the spatial and temporal 
bandwidth constraints relevant to astronomical observations. 


This paper is organized as follows. In Section 2, we first review the quantum theory of thermal radiation. 
We use a full multi-temporal mode analysis to derive the density operator of filtered thermal radiation in the 
limit of coherence time Tc « T, the observation time. This condition is typically satisfied for both radio and 
optical astronomical observations. The details of the calculation appear in the Appendix. In Section 3, we 
briefly review the quantum Cramer-Rao bound before applying the results of Section 2 to explicitly evaluate 
it for filtered thermal radiation. We obtain the result that the right-hand side of eq. m is a fundamental 
quantum limit to the relative sensitivity of any unbiased estimator of the source temperature. We close with 
a discussion of the implications of this result in Section 4. 


2. Quantum-mechanical description of input field 
2.1. Filtered thermal radiation 

Regardless of the frequency range, the electromagnetic field input to an antenna or telescope is described 
by a time-dependent positive-frequency field operator E(t) (in units of -^photons • s"^) given by 


E{t) = 

JO 


Qji/ 6 


-i2T:ut 


dv. 


(4) 


Here, the are un-normalized single-frequenc y annihilation op e rators satisfying ajy,a|,,] = - v') 

which implies that \E{t),E^{t')\ = As in ( Lieu et al. 2015 : ZmuidzinasI 2015 1. we are assuming 

the input field to be of a single polarization and in a single spatial mode to focus on the central issue - the 
additional generalit y does not sub stantially alter the result. If the input field is fr om a thermal source, it 
is a Gaussian field ( Shanirol 20091 Sec. Ill) with the quantum expectation yalues ( Mandel and WoH 1995t 
Zmuidzinasll2003a):- 


(a.) = 0, (5) 

{a^av') = Q ( 6 ) 

= nth{v)5{v-v'), (7) 

where nth{v) is the mean occupation number in a thermal state giyen by the Planck formula 

= -, ( 8 ) 

gkTs — 1 

for Ts the source temperature (in K) and fc, Boltzmann’s constant. In terms of field operators, these relations 
imply 


(A(t)) = 0, 

:= {E{t)E{t')) = 0, 

:= {E\t)E{t')) = dv. 


(9) 
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Here, the functions (t, t') and (t, t') are the phase-sensiti ve and (normal ly ordered) phase-insensitive 
correlation functions of the field respectively (cf. eqns. (61)-(62) of Shapirol ( 2009IH . Note that these functions 
(as well as the mean) depend only upon the time difference r := t' -1, indicating the s tatistical stationarity 
of the field. Using the Gaussian moment-factoring theorem ( Mandel and Wolf IQQSl) . all higher moments 
of the field operators can be expressed in terms of these functions, which therefore constitute a complete 
description of the field. 

The thermal field described above is broadband. In practice, a measurement operates on only a finite 
band of the input field that is determined either by insertion of filters or the response of the measuring 
device. We will accordingly assume that the input is passed through a filter centered at frequency 
and with a flat profile over the band [z/q - Ai//2,i/o + Ai//2]. In a typical radio astronomy measurement, 
e.g., we may have z^o = 1 GHz, and Az/ a few MHz. Further, in radio astronomy, the source temperature 
Ts is typically such that hyolkT s « 1, so that the Rayleigh-Jeans approximation to eq. ([8|) may be used 
( Burke and Graham-Smith 2010ll . As a further simplification, we may assume that the mean occupation 
number of eq. m is approximately flat in the band [zzq - Az//2, zzq + Aiyl2] at the value ng - kTsjhvQ. As a 
result of the above assumptions, for the field at the output of the filter, we replace n^y,{v) in eqns. © and 
(HH) by 


i(z/) = 


no = 

0 


kTs 

hUQ 


if z/ € [z/Q - Az//2, z/q + Az//2] 
otherwise. 


( 12 ) 


In this quasi-monochromatic regime, n{v) is essentially the (dimensionless) power spectral density of the 
field, in that multiplication by huo gives the average power per unit frequency (in W-Hz~^) of the field. 

The held at the output of the bandpass hlter continues to satisfy eqns. (O and (fTUl) but Eqn. m is 
modihed to 


= uq ■ Ah' sine [Azzr] e 


-i27TUoT 


(13) 


wher e sine (a:) = sin(7rx)/(7rx) is the sine function. The coherence time Tc of the output radiation ( Mandel 


1959 ). which is the approximate “width” of |Ar^"^(r)|, is then 


Tc 




A:(")(r) 

is:(")(o) 



(14) 


and is also the separation between the peak of |Ar*'"^(T)| at r = 0 and the hrst zero. For the radio frequency 
example above, we have Tc ~ 1 fjs. In the optical regime, Tc is typically much smaller even downstream of an 
optical hlter, of the order of nanoseconds. 


The measurement on the held takes place in a hnite time interval T of duration T s, which we take to be 
[-T/2,T/2]. In the remainder of the paper, we assume that the parameters z/q, Azz (or Tc), and T are hxed 
and known, while rzg (or equivalently, Tg) is the single unknown parameter that we wish to estimate. We 
further assume that Tc « T, which is usually the case for the detection of faint sources, for which T could 
range from seconds to hours to days. 


2.2. Modal description 

Quantum information and metrology usi ng electromagnetic held s (Helstrom 19761: Holevo 2011 ). and 
Gaussian quantum information in particular (|Weedbrook et al.ll2012t IOlivaresll2012i) . is usually set up and 
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studied on a finite set of modes that are excited in the problem under consideration. In order to obtain 
a modal description of the field from that in terms of its mean and correlation functions in the previous 
subsection, we expand the field in the measurement interval T in terms of a complete orthonormal set of 
positive-frequency “Fourier-series” traveling-wave modes given by 


0m(t) 


^expH^) 

0 


if t € [-T/2,T/2] 
otherwise, 


for TO = 0,1,2,.... The modal annihilation operators 

/ oo 

E{t)(j)ln_{t)dt, TO =0,1,2,... 

CX3 

then satisfy the (normalized) commutation relations 

~ dmn- 


(15) 


(16) 


(17) 


As in standard quantum mechanics, the field state is described by a density operator p (a positive semidefinite 
operator with Tr p = 1) on a Hilbert space H = where Hm is the (infinite-dimensional) Hilbert 

space of the rn-th mode. Our next task is to obtain p. 


Since E {t^ is a Gaussian field, the density operator p is in a so-called Gaussian state ( Shanirol 2009t 
Holevolboilll . In order to define Gaussian states, we need some notation. For each to, we define the 


quadrature operators 


Qm ~ 


am + al 


V2 


"? Pm — 


- at. 


V2i 


(18) 


satisfying the canonical commutation relations [qm,qn] = [Pm,Pn] = 0 and [qm,Pn] = iSmn- Gonsider the 
vector R e (i?i,i? 2 , ■■■)■= {qi,Pi,q 2 ,P 2 , ■ • ■) of quadrature operators. The mean vector R in the state p is 


R:=(R)^ = (Tr(pRi),Tr(pR 2 ),...) 
and the covariance matrix cr has the (i,_))-th matrix element 

~ ~ ^ ~ ~ ^i))p ~ RelTr p (i?i - Ri) {Rj - 7?^)} . 

With this notation, a Gaussian state is a state whose Wigner characteristic function 

XpU) ■= Tr [p exp Cl R)] 
of the Gaussian form ( Olivares! 2012!) 


is 


Xp {£,) = exp {^-i$,^ R - j , 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


where ^ ■ • ■) is the vector of Fourier variables corresponding to the phase-space coor¬ 

dinates of each mode. The matrix 


0 = ^0; 

m 


( 23 ) 
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is block-diagonal in the 2x2 blocks 


with one block per mode. 


ijj = 


C, :)■ 


(24) 


From eq. (HID, we see that a Gaussian state is completely described by its mean and covariance matrix. 
To calculate these for filtered thermal radiation, we must compute the first- and second-order moments of the 
form (om), {amOra), and for each m and n. In the Appendix, we present the detailed calculations 

for the above quantities in the long observation time limit Tc « T. The results are 


(am) “ 0 

{dm dn ) ~ 0 

' jyp _ Hi 

(a^an) - norect 


■ Sr, 


where 


rect(a;) = 


1 if |a:| < 1/2 
0 otherwise 


(25) 

(26) 

(27) 


(28) 


is the rectangle function. In other words, in this « T limit, we have (a])ja„) ^ 0 for m + n and the 
average number of photons 


{alram)-no for m e AI = 7", I,---, 


(29) 


i.e., in M := TAi/ “approximately single-frequency” modes that are within the bandwidth of the filter. 
Together with eqs. dAiD-dMD, eq. (12911 implies that the covariance matrix cr of p is 


© (^*d © (i»). 

nc A/f ' ^ ' rn^Ai 


(30) 


m^Ai 


where I 2 is the 2x2 identity matrix. This implies that M = AvT of the modes (those in the set M. defi ned 
in eq. (l29l) l are in independent thermal states each with no photons on average given by (IOlivaresll2012ll :- 

Pth(?^o) = ^ . E ( ^° 1 ] \k){k\ = f exp(- — \\a){a\(fa. (31) 

’T'O + 1 fc=o ^ *^0 + 1 / 7rno JC \ no/ 

In the first representation, pth(no) is a mixture of number states with a Bose-Einstein distribution, while in 
the second, it is a zero-mean circularly-symmetric Gaussian distribution of coherent states {|q;)}, i.e., of the 
eigenstates of the annihilation operator - a|a) = a |a) - of the relevant mode. The remaining modes are all 
in the vacuum state |0) (0|. Thus, the overall state is 


\meA1 / / 


(32) 


Since the modes m i M carry no information about no, it is sufficient to make measurements on just the 
modes in the set M, effectively reducing the problem to one involving a finite number of modes. 
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3. The quantum limit on estimating source temperature 


3.1. The Quantum Cramer-Rao bound 


The classical Cramer-Rao bound ( Cramei 1946; R.aol 1945 : Van Trees 2001) provides a lower bound on 
the variance Var 9 of any unbiased estimator 9 of an unknown parameter 9 indexing a family of probability 
distributions {Pg} on a given sample space. An unbiased estimator of 9 is one that satisfies E[0] = 9, where the 
statistical expectation value is taken with respect to Pg. The Cramer-Rao bound is widely used to provide 
limit s on and benchmarks for the performance of communication and measur ement systems (Van Tre€^ 
20011) . It has also found application in the design o f astronomical inst ruments (jZmuidzinasI l2003bl) and in 
measurements of the cosmic microwave background (|Yadav et al.l 120071) . 


In a quantum estimation problem, instead of probability distributions {Pg}, we are provided with a 
family {pg} of density operators of a given quantum-mechanical system depending on the unknown parame¬ 
ter. The additional feature of the quantum estimation problem over its classical counterpart is the freedom 
of choosing the quantum measurement that generates a probability distribution from pg. All possible quan¬ 
tum meas urements can be mathematica lly described by an object called a positive-operator-valued measure 
(POVM) ( Helstrom 19761 : Holevo 2011 ). which subsumes the well-known observables of standard quantum 
mechanics. T he quantum Cramer-Rao bound (henceforth “q-CR bound”) ( Helstrom 1967 . 19681 1973 . 19761: 
Holevol 2011 ) provides a lower bound on the variance of any unbiased estimator 0 of a parameter 9 indexing 
a family of density operators {pg} optimized over all possible POVMs subject to the unbiasedness condition. 

W e simply state the result of the q-CR bound here, referring to ( Helstroi:i]ll967 . 19681 19731 19761: Holevol 


20111) for details. We are given a family {pg} of density operators depending on the parameter of interest 9. 
The operator equation 


= 2 i^spg + pgLg ), 


(33) 


has a unique Hermitian solution Lg = Lg when pg has no zero eigenvalues ( Bhatia 2007 ). The operator Lg 
is called the symmetric logarithmic derivative (SLD) in analogy with the classical case. The quantity 


lQ{9)=i:v{pgLl) = {Ll)^ 


is known as the quantum Fisher information and the q-CR bound reads 

1 


Var 9 > 




(34) 


(35) 


and is valid for any unbiased estimator satisfying E[0] = 9, where the expectation is over the probability 
distribution induced by the POVM on the state pg. 


3.2. Estimating source temperature 


The q-CR bound was originally developed in the context of quantum optical communication in the 
years following the invention of the laser, so it was natural for the early work to focus on the experimentally- 
important Gaussian states of light, particularly on the estimation of the mean vector (eqn. (lAll) ') of Gaussian 
states. The pro blem of estimating the average photon number in a thermal state was also considered by 
Helstrom ( 19681) . Very recently, the estimation of a general parameter indexing single-mode and multi-mode 


Gaussian states has been considered by iPinel et al.l (|2013[) and iMonrasI (120131) respectively. 







































































With the density operator for filtered thermal radiati on now at hand from our modal decompositio n 
in Sec. 2.2, we may invoke the required q-CR bound from (lHelstromlll968t iPinel et al.ll2013i lMonrasll2013ll . 
For completeness, however, we re-derive the q-CR bound for estimating the average photon number ng, to 
begin with, in a single-mode thermal state ptuino). From the first representation in eq. (1311) . we see that 


no + 1 \ no + 1 / 


(36) 


where N = o' a = Z^o ^ 1^) (^1 number operator. Since no > 0, pth(no) has no zero eigenvalues and a 

unique SLD exists. To find it, we compute the derivative 


(37) 

(38) 


9pth(no) 

1 

[iV/ no 


dno 

(no-h 1)2 

no \ no -h 1 / 

\no + l) 


= Pth (no) 


N 


1 


Comparing with eq. 


no(no + l) no-hi 

we obtain the SLD operator 

N _ 1 

no(no-hl) no-h 1 


Lno - 


(39) 


Note that L„p commutes with pth(no), reflecting the fact that the {pth(no)} commute with each other. The 
quantum Fisher information is 


Igino) = {Llg) 


Pth(»o) 


1 


(no -h 1)2 
1 

(no -h 1)2 
1 


iV2 

T 


2N 


-h 1 


n, 
2{N) 


0 ’^-o 

2 


Pth(tio). 

Pth(rao ) ("^^PthCno) 


2{N) 


Pth(tio) 


-h 1 


n, 


0 


no 


(40) 

(41) 

(42) 

(43) 


no(no -h 1) ’ 

where we have applied Gaussian moment factoring (Mandel and Wolf 19951) to obtain eq. (1421) - cf. also eqn. 
(25) of ( Helstrom 1968 ) and eqn. (19) of ( Pinel et al. 2013I) for the final result. 

According to eq. (15^ . the input field is the tensor product of M = AnT modes each in the thermal state 
Pth (no). From the additivity of the quantum Fisher information for tensor-product states (which follows 
directly from eqs. (I551) - (IM1) L we get the total Fisher information 

AvT 


rtotal 


(no) = 


leading to the sought q-CR bound 


Var ho > 


no(no -h 1) ’ 

no(no -h 1) 

AvT 


(44) 


(45) 


valid for any unbiased estimator hg for ng. The relative sensitivity of any unbiased estimator hg therefore 
satisfies 


Var ho ^ ng -h 1 


noAvT 


(46) 
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This lower limit on the relative sensitivity of any unbiased estimator of no is our main result. From eq. m, 
we see that the relative sensitivity of an unbiased estimator Tg of the source temperature similarly obeys the 
limit 


Varfg ^ 1 + ^ 

T2 ■ Ai/T 


(47) 


4. Discussion 


First, we note that the estimators corresponding to photon co unting and the est imator 


corresponding to the two-detector scheme of iLieu et all (1201511 are unbiased (IZmuidzinasll2ni,5ll . therefore the 
q-CR bound deriv ed above applies to them. Eq. dUl) agrees with the relative sensitivity of eq. found by 
Zmuidzinas (2015f) for both estimators. It also agrees with the sensitivitv of ea. m for the ideal radiometer in 


the photon-rich hq » 1 regime. This is not surprising because in this limit, the radiometer - which performs 
a heterodyne measurement followed by post-processing to convert inf ormation on the two quadratures into a 
photon-number or energy measurement - essentially counts photons (|Nitvanandalll994l : IZmuidzinasll2003ar) . 
However, the limit of eq. (H51) is greater than that of eq. (I2|) in the Tsamp « « T regime for which the 

benefit of the two-detecto r scheme is claimed. As such, it provides a refutation of eq. ([2]) independent of 
that in ( ZmuidzinasI 2015 ). 


The q-CR limit eqs. (H51) - (IT7l) says much more, however. Recall that the q-CR bound is applicable 
to all POVM measurements made on the input state subject to the unbiasedness condition. Since any 
concrete measurement scheme, ideal or non-ideal, corresponds mathematically to a POVM, we have shown 
that no possible unbiased measurement can improve on the sensitivity of eqs. (|46l)-g2l). This sensitivity is 
therefore “future-proof” and cannot be improved upon by an ingeniously designed measurement or by future 
technological developments. Interestingly, when no » 1, the ideal (noiseless) radiometer already approaches 
this sensitivity limit, making it near quantum-optimal. 
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A. Calculation of mean and covariance matrix 


We calculate here the first and second-order moments {am), {am an), and (a)„ a„) 
radiation described in Section 2.1. Using eqn. (0, we get {am) = (aL) = 0 so that 


for the filtered thermal 


R= 0. 


(Al) 


Similarly, eq. m gives 


{cirn ^n) ~ 0- 


(A2) 





























10 


Using eq. (nsi), we get 

1 rTII ^T/2 , ^ / 

(*’'') “p(-‘ 


27r(mt - nt 


T 




dtdt' 


^ 1-1'“^ /<r!'|, " (^)l) 


no 
2rcT J-3 


-(T-ICI) 
, exp 


j / -o r m+n] ] 

/ aCexp -z— (m-njC / arsine — exp -z27r z/q -n 

J-T 2r^ ^^JJ_(T-ici) \tJ L 2T \ ) 


(A3) 

(A4) 

(AS) 


where we have changed variables to (^ = t' + t and t = t' -t. We now exploit the fact that the coherence time 
Tc « T, the observation time. Consider a time duration ct^ a few coherence times long, where c is a small 
number, say c 5. We split the integral in eq. (ESI) into three parts:- 


(oLan)= 


no 


2tcT J-t 


L 


-(T-cre) 


r 

f 


np 

2TcT J-(T-crc) 
no 


dC, exp 
dC exp 


-z- 


. 27r 


(m-n)cj J 


T-ICI 

(T-ICI) 

T-ICI 


dr sine 


(i)exp(-i2. 


2tt , ,1 r^-KI , . /t\ / „ 

i — (nz-n)C / ar sme — exp -z27r 

2T^ ^'’Jj_(T-|ci) \tc) V 


I/O - ■ 


vp 


m + n 


2T 
m + n 
2T 




2TcT Jr-CTa 


d( exp 


r 27r 1 

cT-ICI / 

T \ / 

r m + nl 

-z—(m-n)C 

/ dr sine 

— exp -izTT 


L 2T J 

J-(T-K\) \ 

Tc) \ 

L 2T J 


.")• 


(A6) 


Consider the first and third terms above. For any the inner integrals in these terms range over an interval 
of size less than or equal to 2cTc. Since, sinc(x) < 1, the inner integral is bounded in absolute value by 2 cTc. 
The outer integral is over a ^-range of cTc, so each of these terms is bounded in absolute value by npc^Tc/T. 
Since we are assuming « T, we can neglect these terms in this limit. On the other hand, the inner integral 
in the second term ranges over a r interval around r = 0 that is equal to or greater than cTc- Since most of 
the area under the sine function is contained in the first few sidelobes, we can approximate, for |C| < T - cTc, 


where 


X 


dr sine 


(T-ICI) 

dr sine 


i„c(i)exp(- 

(i)pxp(-.2 


l/Q- 


m + n 


Vp - ■ 


2T 
m + n 
2T 




/) 


= Tr rect 


Vp ■ 


m+n 

2T 


Al 


rect(x) 


1 if |x| < 1/2 
0 otherwise 


is the rectangle function. Evaluating the outer integral, again using Tc « T, gives 

^0 'J' 


{al^an) - zzorect 


Av 


■ Sr, 


(A7) 


(AS) 


(A9) 
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